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Abstract 

Suggested holographic duals of QCD, based on AdS/CFT duality, predict that one should be 
able to vary the scales of colour confinement and chiral-symmetry breaking independently. Fur- 
thermore they suggest that such independent variation of scales can be achieved by the inclusion of 
extra 4-fermion interactions in QCD. We simulate lattice QCD with such extra 4-fermion terms at 
finite temperatures and show that for strong enough 4-fermion couplings the deconfinement tran- 
sition occurs at a lower temperature than the chiral-symmetry restoration transition. Moreover 
the separation of these transitions depends on the size of the 4-fermion coupling, confirming the 
predictions from the proposed holographic dual of QCD. 
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I. INTRODUCTION 



The AdS / CFT duality of Maldacena , which indicates that Af = 4 supersymmetric 
Yang-Mills field theory in 3+1 dimensions is dual to M = 8 super-gravity in 4+1 dimensions, 
inspired people to search for a gravity/string theory which is dual to QCD. Since such duality 
would map the strong-coupling regime of QCD to the weak-coupling regime of a classical 
gravity theory, this promises to greatly simplify certain QCD calculations. For a review of 
the early work in trying to construct such holographic duals for QCD without quarks, we 
refer the reader to the review article by Aharony [3]. The inclusion of quarks was addressed 
by Karch and Katz 41. Another such model for including quark flavours was suggested by 

(| and studied further by Antonyan, Harvey, Jensen and Kutasov 

These papers 



]) observe that 



Sakai and Sugimoto 

j?! and Aharony, Sonnenschein and Yankielowicz 
the scales of confinement and chiral-symmetry breaking can be varied independently in 
these proposed QCD duals, and that the chiral-symmetry breaking scale must always be 
shorter than or equal to the confinement scale. Note that all these attempts at constructing 
holographic duals of QCD only claim validity for large N c . 

Measurement of the scales of confinement and chiral symmetry breaking can be addressed 
in lattice QCD by measuring the deconfinement and chiral-symmetry restoration temper- 
atures for hot QCD. For quarks in the fundamental representation of SU(3) colour, these 
two transitions have been observed to be coincident, within the limitations of lattice mea- 
surements {9]. For this reason there have been a number of attempts to explain why this 



should be so 1(3 . 



Others have argued that the scales of confinement and chiral- 



symmetry breaking could easily be different 



14J . However, it has been observed that for 



quarks in the adjoint representation [15j|, and probably for quarks in the sextet representation 
of colour [16j, the deconfinement temperature is significantly lower than the chiral-symmetry 
restoration temperature. This suggests that the real reason why these transitions appear 
coincident for fundamental quarks, is that the interaction between fundamental quarks and 
antiquarks is too weak to produce a chiral condensate at distances shorter than the confine- 
ment scale. The condensate is then produced at the confinement scale since confinement 
requires chiral-symmetry breaking [lT], [Q . This would then make the two scales and hence 
the two transition temperatures identical. 

The work of Antonyan, Harvey, Jensen and Kutasov [7|] suggests that the introduction 
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of extra 4-fermion interactions in QCD would allow the separation of the confinement and 



chiral- symmetry breaking scales. This ma 



Gross-Neveu 



£es sense, since it is known that models of the 



19] /Nambu- Jona-Lasinio 



20 



21] type exhibit chiral-symmetry breaking with- 



out confinement for strong enough coupling. Hence one might expect that the addition of a 
large enough 4-fermion coupling to QCD could produce chiral-symmetry breaking without 
confinement. We therefore choose to study lattice QCD with extra 4-fermion interactions 
of the Gross- Neveu/Nambu- Jona-Lasinio type at finite temperatures, looking for evidence 
for separate deconfinement and chiral-symmetry restoration transitions. Here we have been 
careful to make sure that the 4-fermion coupling is not so strong as to produce sponta- 
neous chiral-symmetry breaking without gauge fields, since this phase is separated from the 
chiral-symmetry restored phase by a bulk transition. 

Early studies of such a model showed some evidence for separate transitions for strong 
enough 4-fermion interactions [22], but the lattice size (8 3 x 4) was so small that the ef- 
fect could have been a finite lattice size artifact. We study this model for large 4-fermion 
couplings on 16 3 x 4 and 24 3 x 4 lattices, and for intermediate 4-fermion couplings, on 
2 2 x 24 x 4, 24 3 x 4 and 32 3 x 4 lattices (the 12 2 x 24 x 4 'data' is from a previous study 



23l|) . For strong 4-fermion couplings, we find that the deconfinement transition takes place 



at a much lower temperature than the chiral-symmetry restoration phase transition. At 
intermediate 4-fermion couplings, the 2 transitions are much closer, but still clearly sepa- 
rate. At weak 4-fermion couplings, previous work has indicated that the two transitions are 
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23J. 



coincident within the limits of our simulations 

In section 2 we introduce the action for lattice QCD with additional 4-fermion interactions 
and discuss some of its properties, and how it is simulated. Section 3 presents our simulations 
at zero gauge coupling which are necessary to identify the limits on the 4-fermion couplings 
we can use. In section 4 we describe our simulations and present results. Section 5 discusses 
our results, draws conclusions and indicates directions for future research. 



II. X QCD 



The lattice QCD action with additional 4-fermion interactions which we choose is one 
we have previously called xQCD 22] . This is based on the continuum Euclidean space-time 
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Lagrangian density 

C = \f, u F^ + 4,{If)+ rr# - ^[(^) 2 - fcrs^) 2 ]. (1) 

Introducing auxiliary fields a and 7r yields a new equivalent Lagrangian density which is 
quadratic in the fermion fields. 



C = -F^F^ + ij(p+ a + z7T7 5 r 3 + m)il) + T^V + A (2) 



The chosen 4-fermion term is of the Gross-Neveu/Nambu-Jona-Lasinio type. It breaks 
flavour symmetry, but preserves the reduced chiral symmetry of the staggered-fermion im- 
plementation of quarks on the lattice. Unlike the action which includes Vi' and 'r 2 ' terms 
and preserves the full chiral flavour symmetry, it has a real positive fermion determinant 
which is essential for lattice simulations. The terms involving the chiral auxiliary fields (a, 
7i") remain local in the staggered-fermion lattice transcription. 

The lattice version of the action is 

S = 0£ f 1 " ^M^nUUUU)} + £ \N fl (o 2 + 7T 2 ) 

n L O J - O 



Nf/A 

+ E Ex/ 

/=i s 



ft + m + 77 E(° r * + ^ e7r *) 
lb 



X/ (3) 



where ^ is the standard staggered gauge-covariant ft, e = (— \Y+v+ z + t ) s ruri s over the sites 
of the lattice, s runs over the sites of the dual lattice on which the auxiliary fields reside 
and % runs over those sites of the dual lattice adjacent to the site of the lattice on which the 
fermion field resides. This transcription preserves the U(l) chiral symmetry of the staggered 
fermion formulation. 7 = 12/ A 2 [31]]. 



We simulate this action using the RHMC algorithm 24j], where the fractional powers of 
the fermion determinant, required when the number of flavours Nf is not a multiple of 8, 
are obtained using a rational approximation to the fractional powers of the quadratic Dirac 
operator, and global Metropolis accept/reject steps make the algorithm exact. For subtleties 
associated with applying the RHMC algorithm to this action, see reference j^. As noted 
in our earlier work, the presence of the 4-fermion interaction makes the Dirac operator non- 
singular in the chiral m = limit. We make use of this fact to simulate at m — 0, where 
there is an exact £7(1) chiral symmetry, so that the chiral-symmetry restoration occurs at a 
true phase transition at and beyond which the chiral condensate vanishes. There is, however, 



no true order parameter known for the deconfmement transition in QCD with dynamical 
fundamental quarks, which is seen as an abrupt jump in the Wilson Line (Polyakov Loop). 
Hence this transition need not necessarily be a true phase transition. 

When we simulate at m — 0, the direction in which the U(l) chiral symmetry is broken 
is arbitrary. The chiral condensate is some linear combination of (ipip) an d ^(V'Ts^sV') or 
of a and 7r (£5, the flavour version of 75 is the equivalent of T3 for the 4 flavours described 
by a single staggered quark field). On a finite lattice, the direction defined by this chiral 
condensate rotates during the simulation, and the condensates average to zero since there is 
no spontaneous symmetry breaking in a finite volume. We therefore define an approximate 
order parameter, which approaches the true order parameter when the volume of the lattice 
becomes infinite. Experience is that \J1p1p 2 — ^Isi^ 2 or \J o 2 + it 2 are good choices. Here 
the quantities ipifj, etc. are lattice averages for a given gauge and auxiliary field configuration. 
These choices lead to some finite size rounding of the phase transition on finite lattices. From 
now on we shall denote these two versions of the chiral condensate as u {ipip)" and "(&)"■ 

All our simulations have been performed with two flavours of massless quarks (Nf = 2). 

III. GROSS-NEVEU/NAMBU-JONA-LASINIO MODEL SIMULATIONS 

To determine the relevant values of 7, we perform simulations with the pure Gross- 
Neveu/Nambu-Jona-Lasinio model, i.e. the theory described in the previous section, only 
without gauge fields. For these simulations, we run at zero temperature. Note that Nf = 2 
QCD with extra 4-fermion interactions reduces to a 6-flavour Gross-Neveu/Nambu-Jona- 
Lasinio model (2 flavours x 3 colours) when the gauge coupling is set to zero. We are 
interested in determining the 7 value for the bulk transition at which the chiral condensate 
vanishes. Since, for our simulations with gauge fields, we want to determine the two finite 
temperature transitions, we need to simulate with 4-fermion couplings weaker than that at 
the bulk transition, i.e. for 7 > 7 C where 7 C is 7 at the bulk chiral-symmetry restoration 
phase transition. In our earlier work, we simulated this theory with only 4-fermion couplings, 
on an 8 4 lattice. For this project we use a 12 4 lattice. 

We performed runs for 1.0 < 7 < 2.5. For each 7 we generated 500,000 length 2.5 
trajectories. Figure [T] shows the chiral condensates (ipip) an d (&) as functions of 7 over this 
range. From this we estimate 7 C ~ 1.7. Comparing this graph with that for an 8 4 lattice 
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from our previous work, we observe that, as we approach 7 C from below, the finite size 
effects only become appreciable when we get very close to the transition, as expected. The 
behaviour of (a) can be understood from that of (ipip), since in the limit of infinite volume, 
these are related by 



This equation remains true when we reintroduce the gauge fields. 

IV. FINITE TEMPERATURE SIMULATIONS AND RESULTS. 

We now turn to simulations with finite gauge couplings. The quark mass m = for 
all these simulations, so that the chiral transition is a true phase transition. Since the 
separation of scales of confinement and chiral-symmetry breaking is expected to occur for 
strong 4-fermion interactions, we start with a large 4-fermion coupling. We choose 7 = 
2.5 > 7 C ~ 1.7. While this represents a large 4-fermion interaction it is not too close to the 
bulk transition when the gauge coupling is taken to zero. Ideally one would like to keep 
the gauge coupling fixed so that 7 is fixed in physical as well as lattice units, changing the 
temperature by changing the temporal extent (in lattice units), N t , of the lattice. As is 
always the case (unless one uses anisotropic lattices) this is impractical, and we fix N t and 
vary the temperature by varying the lattice spacing in physical units by varying (3 = Q/g 2 
where g is the gauge coupling constant. Since we keep 7 fixed in lattice units, 7 changes in 
physical units as we vary (3. For these simulations we use N t = 4. 

Earlier simulations on 8 3 x 4 lattices at 7 = 2.5 suggested that there could be two 
separate transitions. However, on such small lattices, finite size effects are large enough that 
definite conclusions are suspect. We have now simulated at 7 = 2.5 on 16 3 x 4 and 24 3 x 4 
lattices. At each j3 we ran for 50,000 length 1 trajectories, except for (3 = 5.545, close to 
the deconfinement transition, where we ran for 100,000 length 1 trajectories. In figure [2] 
we plot the Wilson Line (Polyakov Loop), which signals the deconfinement transition and 
the chiral condensate, which is the order parameter for chiral-symmetry breaking, against 
(3. The deconfinement transition, marked by a very rapid increase of the Wilson line from 
near-zero, occurs at a much smaller (3 than the chiral phase transition, where the chiral 
condensate vanishes and chiral symmetry is restored. In the region of the deconfinement 
transition, we use Ferrenberg-Swendsen reweighting [26] to determine the position of this 
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FIG. 1: The chiral condensates {ipifj), (cr) for the 6- flavour Gross-Neveu/Nambu- Jona-Lasinio (zero 
gauge coupling) limit of QCD with extra 4-fermion couplings, as functions of 7. 

transition from the Wilson line susceptibility peaks. Examining the distribution of plaquette 
values for each /3 in the range 5.50 < (3 < 5.57 suggests that the deconfinement f3, (3d ~ 5.545. 
Ferrenberg-Swendsen reweighting from (3 = 5.545 yields (3d = 5.547(3). The chiral transition 
is considerably less well measured, since there are much larger gaps between consecutive /?s 
in its neighbourhood, and the finite size effects are large. Examining the 'data' from the two 



lattice sizes and trying to take into account the finite size effects leads us to (3 X = 7.0(2) as 
our estimate for the (3 at the chiral phase transition. Thus, for 7 = 2.5 the deconfinement 
and chiral-symmetry restoration temperatures are far apart. 
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FIG. 2: Wilson line and chiral condensate as functions of j3 for 7 = 2.5 in lattice units. 



We also performed simulations at an intermediate 4-fermion coupling 7 = 5. Here, earlier 
simulations on 8 3 x 4 and 12 2 x 24 x 4 lattices had failed to indicate whether there was one 
transition or two. We have performed simulations on 24 3 x 4 lattices. To clarify chiral 
symmetry restoration, we also performed simulations on 32 3 x 4 lattices close to the chiral- 
symmetry restoration transition. For the 24 3 x 4 lattice we have run for 100,000 length 1 
trajectories for each (3 in the range 5.415 < (3 < 5.460 and for 50,000 for those /3s outside 
this range. On the 32 3 x 4 lattice we ran for 100,000 length 1 trajectories for each of three 
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/3 values. Figure [3] shows the chiral condensate and Wilson line from these new simulations 
along with the old results on 12 2 x 24 x 4 lattices. A gain we find evidence for two separate 
transitions. Close to the deconfinement transition, we find that the finite size effects are 
very small, as was observed at 7 = 2.5. This is reassuring, especially in light of the fact that 
our old 12 2 x 24 x 4 simulations used the inexact R algorithm, while the new simulations 
used the exact RHMC algorithm. We estimate that the deconfinement (3, (3d = 5.420(4). 
Our estimated (3 for the chiral transition is j3 x = 5.450(5). 
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FIG. 3: Wilson line and chiral condensate as functions of (3 for 7 = 5 in lattice units. 



We know from our earlier work that if 7 is increased to 7 = 10, the deconfinement and 
chiral transitions appear to be coincident. The closeness of [3d and (3 X at 7 = 5 compared 
with 7 = 2.5 suggests that the two transitions come together at a 7, a little above 5. 
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V. DISCUSSION AND CONCLUSIONS 



We have simulated lattice QCD with extra 4-fermion interactions of the Gross- 
Neveu/Nambu-Jona-Lasinio type at finite temperatures, and have observed that, for strong 
enough 4-fermion couplings, the deconfinement and chiral-symmetry restoration temper- 
atures are different. Increasing the 4-fermion coupling increases the separation of these 
two transitions, and the deconfinement temperature is always less than or equal to the 
chiral-symmetry restoration temperature. This agrees with the predictions from proposed 
holographic duals of QCD [7, 8]. 

These results can be easily understood, since Gross-Neveu/Nambu-Jona-Lasinio models 
without gauge fields can exhibit spontaneous chiral-symmetry breaking, but do not confine 
the fermions. It also indicates that the reason why for QCD without such terms, the two 
transitions appear coincident is merely that the eigenvalue of the quadratic Casimir operator 
for the fundamental representation of SU(3) co i our is too small for the QCD interactions to 
produce chiral-symmetry breaking at distances less than the confinement scale. Confinement 
then forces chiral symmetry breaking and the two scales are identical. 

When the two scales are different, chiral symmetry is broken on both sides of the de- 
confinement transition. Hence, at this transition, the gauge fields see 'constituent quarks', 
i.e. quarks with dynamical masses produced by spontaneous chiral-symmetry breaking, 
rather than the massless 'current quarks'. These are less effective at screening colour than 
the 'current quarks' and the deconfinement temperature increases. We expect that in the 
limit of large 4-fermion couplings the dynamical quark masses will increase without bound 
and the deconfinement temperature will approach its quenched value. This explains why 
the transition coupling (3d and hence temperature Td vary over a limited range as the 4- 
fermion coupling is varied. These arguments suggest 5.25 <■ (3d ?S 5.6925, where the lower 



limit is from simulations without the 4-fermion interaction 



281 ] and the upper limit is 



the quenched value 29j, |30|], over the whole range of 7 (the inverse 4-fermion coupling), 
7 C < 7 < 00, where j c « 1.7 is the bulk transition 7 when the gauge couplings are switched 
off. On the other hand, the chiral-symmetry restoration temperature T x and coupling (3 X 
will vary from T x = Td and (3 X = (3d for 7 > 70 (weak 4-fermioncoupling) to 00 as gamma is 
decreased to 7 = j c , when chiral-symmetry is always broken [32j. Our simulations indicate 
that 70, the 7 value for which the two transitions coalesce, lies between 5 and 10, and is 
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probably closer to 5. 

We have simulated at 7 = 2.5 (strong 4-fermion coupling) and 7 = 5 (intermediate 4- 
fermion coupling) and previously at 7 = 10 and 7 = 20 (both weak 4-fermion couplings) 231 ] . 
The deconfmement /3s are /3 d (2.5) = 5.547(3), (3 d (h) = 5.420(4), ^(10) = 5.327(2) and 
/3d(20) = 5.289(1), consistent with the above bounds, and increasing as expected with 
increasing 4-fermion coupling. The (3s for the chiral transition are /3 X (2.5) = 7.0(2), 
P x (5) = 5.450(5), /3 X (10) = (3 d (10) and /3 X (20) = p d (20). For 7 = 2.5 (strong 4-fermion 
coupling) crude estimates based on quenched running of the coupling constant, from earlier 
lattice simulations, indicate that T x is an order of magnitude larger than Tj, while for inter- 
mediate 4-fermion coupling, 7 = 5, 2-flavour 2-loop running of the coupling constant yields 
T x « 1.041V 

In our model, the 4-fermion interaction is completely local and thus irrelevant in the 
renormalization-group sense. Thus it will vanish in the continuum limit when the lattice 
spacing goes to zero, and we will no longer have two separate transitions. On the other 
hand, the 4-fermion interaction implied by the proposed holographic dual to QCD is non- 
local, which softens ultra-violet divergences introducing the possibility that it might have 
a non-trivial continuum limit. This point needs further investigation. Even if this does 
not happen, and QCD with extra 4-fermion interactions is only defined with an ultra-violet 
regulator such as the one provided by the lattice, it is a useful model since it allows one to 
study confinement and chiral-symmetry breaking independently in this ultra-violet regulated 
(effective) theory. 

We have restricted ourselves to N t = 4 for this preliminary study. A more complete 
study would require extending this to larger N t . At N t = 4, we would need more (3 values to 
accurately pinpoint the chiral phase transition and to determine its nature. Additional work 
would be needed to understand the deconfmement transition. The fact that the vicinity of 
the deconfmement transition shows little finite-size dependence, makes it likely that, if it is 
indeed a phase transition, then it is first-order. 

A lattice analysis of the non-local Nambu-Jona-Lasinio model suggested by the proposed 
string/gravity dual to QCD should be considered. It is a 4-dimensional non-gauge the- 
ory which exhibits chiral-symmetry breaking but not confinement. If it has a non-trivial 
continuum limit, it represents a new class of 4-dimensional field theories. The next step 
would be to include it in the lattice QCD action, just as we have done with the local Gross- 
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Neveu/Nambu-Jona-Lasinio model, or alternatively to study the 5- dimensional gauge theory 
with right- and left-handed quarks pinned to separate 4-dimensional branes, which produced 
it. 

The methods of lattice gauge theory and extensions of AdS/CFT duality to QCD and 
similar quantum field theories can be used to complement one another in the understanding 
of such theories. The particular example described in this paper shows how these ideas can 
be applied to further the understanding of confinement and chiral-symmetry breaking in 
QCD. 
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